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Abstract 

We study localized modes on a single magnetic impm'ity positioned in the bulk or at the surface 
of a one-dimensional chain, in the presence of a magnetic field B acting at the impurity site. The 
strong on-site nonlinear interaction U between two electrons of opposite spin at the impurity site, 
modelled here as a nonlinear local term, and the presence of the external field induce a strong 
correlation between parallel and anti-parallel spin bound states. We find that, for an impurity in 
the bulk, a localized vector mode (with up and down spin components) is always possible for any 
given value of U and while for a surface impurity, a minimum value of both, U and B is needed 
to create a vector mode. In this case, up to two localized modes are possible, but only one of 
them is stable. The presence of the surface seems to destabilize the bulk mode in the parameter 
region U ^ B, creating a "forbidden strip" region in parameter space, bounded hy U = B + V and 
U = B — V , approximately. 

PACS numbers: 71.55.-1, 72.10.Fk, 73.40. Gk 
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I. INTRODUCTION 



Nonlinear effects are tfie subject of intensive, ongoing researcli in several fields that include 



„>olacu,a. c.y=t*Q,a, Jo.e.h=on-iu„ct.o„ an.y=H, Fe_etic ■natenal.M, photon, 
crystals [5|], photonic lattices [6|], Bose-Einstein condensates in magneto-optical traps [7|] and 



nonlinear metamaterialsjS!], to name a few. This interest is due in no small part to the wide 
range of potential applications to the design and operation of optoelectronic devices. 

In condensed matter, nonlinear effects arise from at least two different sources. One 
of them is a possible strong coupling between an excitation and local vibrational modes. 
In the approximation where one assumes a rapid readjustments of the local vibrations to 
the presence of the excitation, one quickly arrives at some version of the discrete nonlin- 
ear Schrodinger (DNLS) equation, whose exact form depends on the type of anharmonicity 
(or lack of it) of the underlying oscillators 9(] . When the strong interaction tales place 
only at some few sites, or at a single impurity site, it is possible to compute the forma- 
tion of bound states in closed form in one-dimensional [l^, two- dimensional j 1 1 1 and three- 



dimensional lattices [12|, by using a direct extension of the lattice Green function formalism. 

The other source for nonlinearity comes from electron-electron interactions in nanoscale 
devices, such as quantum dots and few impurity models Roughly speaking, the Coulomb 
interaction gives rise to a nonlinear term in the Schrodinger equation, that is usually mod- 
elled as a cubic, nonlocal term in the fermionic field operators. This is usually followed 



by the Hartree-Fock approximation for t 
tive approach in the Coulomb interaction 



le nonlinear term[l^, or the use of a perturba- 



15|. Another recent approach calls for mode^ 



ling 

3- 



the electron-electron interaction by a nonlinear local term in the Schrodinger equation 
While being a sort of oversimplification of the many-body problem, it has the advantage 
of retaining the main features of Coulomb interaction-induced nonlinearities, while allowing 
for speedy computation of quantities of interest for electronic transport. This approach 
has been recently used for a simple computation of the zero-voltage conductance across a 
magnetic impurity [l^. 

In this Letter, we focus on the possible localized impurity modes that can reside on top of 
a nonlinear magnetic impurity, where the source of the nonlinearity is due to strong electron- 
electron interaction effects at the impurity site. It is important to ascertain the conditions 
under which such bound states exist, since they could, for instance, scatter other extended 
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excitations in the system. The precise control of this scattering could be of importance in 
the control of the transport properties across nanoscale impurity regions. 



II. MODEL 

Let us start by considering the problem of a one- dimensional discrete system consisting 
of two linear chains (leads) joined by a strongly nonlinear magnetic region, where electron- 
electron effects are important. In the interacting region, the electron-electron repulsion is 
modelled by UpoaPo-a, where poo- = iV'o.o-P, where ipo^a is the probability amplitude of finding 
an electron of spin cr =1, | on site n = 0. 

The coupled evolution equations for the probability amplitudes are 



dt 



+ \/(^„+i,l+V^„_i,l)+(5„,o(eo,i + f^l^n,TnV^n,i = (1) 



dt 

where V is the nearest-neighbor coupling parameter, U is the Coulomb repulsion energy 
and eo,t = —B, eg,; = -B is the Zeeman energy shift due to the magnetic field B, which is 
assumed to be appreciable only in the immediate vicinity of the impurity site. 

We are interested in stationary solutions of the type '>pn,ait) = exp(i/?t) ipn,a- This leads 
to the system of nonlinear equations 

At this point an general observation is in order. If we were to consider the impurity site 
also coupled to a local, fast elastic vibrational degree of freedom, then in the limit when the 
local vibration is completely enslaved to the electronic motion, there would be an additional 
term of the form '~i\ipn,u\'^i^n,a in Eqs.([T]), ([2]). In that case, the system would be formally 
equivalent to a birefringent, nonlinear optical Kerr impurity embedded in a weakly-coupled, 
linear waveguide array, in the absence of four-wave effects, and ipn,] and tljn,i would represent 



the electric field amplitudes for the TM and TE mode, respectively [l7l]. 

In order to keep the magnetic and electronic repulsion effects well separated from possible 
polaronic effects, in this work we restrict ourselves to Eqs. (111), (I2l). The interplay of external 
magnetic field, electron-electron repulsion and electron-phonon interaction effects on a single 
impurity, will be described elsewhere. 
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FIG. 1: (Color online). Left: B = case. Amplitude of localized state for both spins as a function 
of nonlinearity. Right: U = case. Amplitude of localized state for both spins, as a function of 
magnetic field. 

III. IMPURITY IN BULK 

Here the chain occupies the interval — oo < n < oo, with the impurity site at = 0. We 
look for a localized mode centered on the impurity site, = A^- ^i"', with < |^| < 1. 
This ansatz leads to the equations: 



On the other hand, from probability conservation, we have 1 = J2n \'^n,^\'^ = ^„ \4'n,i 
This implies the additional equations 



-P + 2Vi^ + + UA\ = 



(3) 



and 




(4) 




(5) 



and 




(6) 



From Eqs.([3]),(jl]),([5]) and ([6]), we obtain two coupled equations for S^a'- 
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FIG. 2: Amplitude of localized mode as a function of nonlinearity and external magnetic field, for 
spin up (left) and spin down (right). Darkest (whitest) shade corresponds to = —1(1). 



and 



B U 



(8) 



It is clear from these equations that, the transformation ^ C,i, B —B leave the 
equations invariant. Since U > due to the repulsive nature of the Coulomb interaction, 
we only need to deal with, say, -B > 0. The behavior for the opposite sign of B is obtained 
by simply exchanging and 

Before analyzing the general case, let us discuss briefly a couple of important special 
cases. 

(a) B=0. In this case, there is no physical distinction between and and Eqs.([7]) 
and ([8]) collapse into a single one: = 1 + {U/V)C,{{1 — ^^)/(l + ^^)), with solution 



U 
2V 



U 

w 



-1, 



(9) 



provided U /V > 2. Since ^ < 0, the spatial mode profile is staggered. Figure 1 shows ^ in 
terms of U/V. As expected, an increase of nonlinearity reduces the width of the localized 
mode. 

(b) U=0. In the absence of Coulomb repulsion, the only source for localization is given 
by the presence of the external field B. The equations read now ^| = 1 — {B/V)C,i,C,J = 
1 + {B /V)^i, with solutions 



B \ 



2V J 



-) +1 

2V J 



(10) 
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FIG. 3: (Color online). Spatial mode profiles for U jV = 3 and varying values of B/V. Top left: 
B/V = 1. Top right: B/V = 2. Bottom left: B/V = 4. Bottom right: B/V = 6. Solid (dashed) 
curves denote the spin up (down) mode. 

and 




Figure 1 shows the allowed as a function of magnetic field. In this case, no minimum 
field strength is needed to create a localized mode, but the mode profile corresponding to a 
spin parallel (antiparallel) to the external field is staggered (unstaggered). 

For the general case, we solve Eqs.([7]) and ([8]) numerically, selecting the real roots that 
lie in the interval (—1, 1). Results are shown on Fig. 2 in the form of a nonlinearity-magnetic 
field phase space diagram, showing the amphtude of for a given value of U/V and 

B/V . For spin up, the amplitude is always negative and increasing with increasing B, 
provided U /V < 2. When U/V > 2, there is a sharp boundary in U — B space separating 
positive amplitude values, in the "low" B/V region, from negative amplitude values, in the 
"high" B/V value. Upon crossing the boundary, in the direction of increasing B, the sign 
of the amplitude jumps discontinuously from +1 to —1. The shape of this critical boundary 
can be computed exactly from Eqs.(l7]) and (IE]): We set = 1 in Eq.(IH]), and solve for 

= -{B/2V) + ^{B/2Vy + 1. Next, we insert this into Eq.d?]), and solve for U in terms 
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of B, obtaining: 

(B\ _ A + {B/Vy-{B/V)y/iB/Vy + A 
\yJc -iB/V) + ^iB/Vy + 4 

For spin down, the amplitude of the locahzed modes is always positive. For U/V < 2, 
the amplitude decreases monotonically towards zero with an increase in B. For U/V > 2 
however, the amplitude starts at some positive value at i? = 0, then increases further with 
an increase in B, until B /V reaches the critical boundary computed above. After that, an 
increase in B leads to an decrement of the amplitude towards zero. 

In Fig. 3 we show some spatial profiles ipn,a for a fixed value of nonlinearity, and increasing 
values of magnetic field. Note that for B > 0, the localized mode for spin down is always 
unstaggered, while for spin up, the mode is unstaggered at first, and becomes staggered after 
certain value of magnetic field, given by Eq. f|T2|) . 



IV. SURFACE IMPURITY 



We now consider the case when the magnetic impurity is at the very beginning of a semi- 
infinite lattice. We relabel the previous chain, so that the first site is now at uq = 0. The 
stationary-state equations read now 

- /3 ^0,. + V^l,^ + ieo,a + U\^o,a\^)^0,a = 0, (13) 
- f3 ^n,a + V{llJn+l,a + ^n-l,a) =0, n = 0, 1, 2, . . . (14) 

We proceed as before and pose a solution of the form = Ao-^", where < |^| < 1 



and n = 0,1,2,.... After replacing this ansatz into Eq. fll3l) and f[Hj) . one obtains (3 = 
V^a + eo,(7 + UAi^ and P = + (V^a), which implies 

= {ejV) + {U/V)Al^- ^^^^ 

On the other hand, from the normalization condition, 1 = I'^n.tP = Yin I'^n-iP' "^^ obtain 
A| = 1 — and Aj = 1 — After inserting this into Eq. (fT5|) and after using eo,cr = —crB, 
one arrives at two coupled transcendental equations for 

^ -iB/V) + iU/V){l~ei)' ^^^^ 
= (B/V) + {U/V) (1 - ■ ^^^^ 




FIG. 4: Left: B = case. Amplitude of localized state for both spins as a function of nonlinearity. 
The solid (dashed) curve denotes the stable (unstable) mode. Right: U = case. Amplitude of 
localized state for both spins, as a function of magnetic field. 

As was done for the impurity in the bulk, let us consider first two special cases. 

(a) B=0. In this case, = ^ ^, where ^ satisfies ^(1 - ^'^) = {U/V)-^. This 
implies that a critical value (U/V)c = {3/2)^/3 ~ 2.6 exists, such that for (U/V) < {U/V)c, 
no localized state exists. For (U/V) = {U/V)c, there is exactly one localized state, and 
{U/V) > {U/V)c, two localized states are possible; for one of them, ^ increases towards unity 
as nonlinearity is increased (unstable mode), while the other decreases ^ as nonlinearity is 
increased (stable mode). Since C/ > 0, this imphes that ^ is also positive, and thus, the 
localized mode is unstaggered. Figure 4 shows ^ in terms of U/V. 

(b) U=0. In the absence of Coulomb repulsion, we immediately obtain = —1/{B/V) 
and = 1/{B/V), where > 1 to ensure < 1. Figure 4 shows the allowed Ct^Cj, ^ 
a function of magnetic field. The mode profile corresponding to a spin parallel (antiparallel) 
to the external field is staggered (unstaggered) . 

Although Figs. 1 and 4 look qualitatively similar, we note some interesting differences: 
For the surface impurity, and in the absence of magnetic field, the amount of nonlinearity 
needed to effect a localized state is higher that in the bulk impurity case. Also, in the linear 
case {U = 0), one needs higher values of magnetic field to create a surface localized mode. 
The surface, or boundary of the system is acting in a repulsive manner. 

For the general case, there are two important regimes, B/V <1 and B/V > 1. 

B/V < 1: In this case, there is a critical nonlinearity value {U/V)c such that, for < 



8 



{U /V) < (f//V^)c, no localized state exists, while for {U /V) > {U/V)c, there are two localized 
states. One of them (stable mode) becomes narrower with increasing nonlinearity; the other 
(unstable mode) becomes wider. 

B /V > 1: In this case, we have three critical values for {U /V), that create four regions: 
(i) < (U/V) < {U/V)ci'. This is the region of "small" nonlinearity, where magnetic field 
dominates and there is always a localized mode, with < 0,^j > 0. As (U/V) {U/V)ci, 



(iii) {U/V)c2 < (U/V) < {U/V)c3: Two localized modes are possible here. One of them 
(stable) becomes narrower with increasing U /V, while the other (unstable) becomes wider 
with increasing U/V. 

(iv) (U/V) > {U/V)c3: The previous unstable state disappears, and only stable mode re- 
mains. As (U/V) oo, both and approach zero. 

Figure 5 shows all these regions in the form of a nonlinear ity-magnetic field phase space 
diagram. The number inside each region denotes the number of localized surface modes (each 
one with a spin up and spin down component). We note, in particular, the existence of a 
region, roughly delimited by the straight lines (U/V) = (B/V) + 1 and (U/V) = (B/V) — 1, 
where no localized states exists. These are obtained as limiting cases of the curves that 
separate regions "1" (top) and "0", and regions "0" and "1" (bottom). For the first case 
the curve is U/V = [{B /V) + !]/[! — (V/B)'^], while for the second case, the curve is given 
by U/V = [{B /V) — !]/[! — {V/BY]. In Fig. 6 we show examples of spatial profiles ipn,a- We 
nota that their staggered/unstaggered character is rather similar to the "bulk" case (Fig.3). 

V. SELFTRAPPING DYNAMICS 

We consider now the dynamical excitation of a localized mode on top of the magnetic 
impurity, comparing the bulk and surface cases. We place initially both electrons (with 
opposite spins) on the impurity site and follow the evolution of their probability densi- 
ties, according to Eq.([T]). In particular, we focus on the time-averaged selftrapped fraction 
remaining the initial site, defined by 



(ii) {U/V),i < (U/V) < {U/V),2- No localized mode exists. 




(18) 
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FIG. 6: (Color online). Examples of surface localized modes. Solid (dashed) lines refer to spin up 
(down). Top left: B /V = l,U /V = 4 (stable mode). Top right: B /V = l,U/V = 4 (unstable 
mode). Bottom left: B/V = 3, U/V = 5. Bottom right: B /V = 3, U/V = 1. 

and examine how {P„) depends on electron interaction U and magnetic field B. 

Figures 7 and 8 show the selftrapping results for the "bulk" and surface cases. While the 
behavior in both cases is qualitatively similar, we note that in general, it is a bit easier to 
selftrap in the bulk case than in the surface case. On the other hand, for the surface case, 
there are sharper boundaries separating the untrapped from the selftrapped regime. Also, 
in both cases and for positive B, selftrapping for the spin up component, is largely inhibited 
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in a wide strip around the region U ^ B, while the spin down component is dominated 
by magnetic field effects. Of course, for [/ « B, nonlinearity effects are not important 
and self trapping increases gradually with magnetic field, as one expects for a linear site 
impurity. On the contrary, for f/ » B, selftrapping is abrupt and its threshold is higher 
for the surface case than for the bulk case. Again, this is a manifestation of the repulsive 
nature of the system surface in one-dimension. 
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FIG. 7: Magnetic impurity in "bulk": Left (right): Time-averaged selftrapped fraction of spin 
up (down) remaining on the initial site, as a function of U and B. Dark (white) shade denote 
selftrapped fraction close to zero (one). 




FIG. 8: Magnetic impurity on surface: Left (right): Time-averaged selftrapped fraction of spin 
up (down) remaining on the initial site, as a function of U and B. Dark (white) shade denote 
selftrapped fraction close to zero (one). 
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VI. CONCLUSIONS 



We have examined theoretically the formation of bound spin states on top of a narrow 
magnetic impurity region in the presence of an external magnetic field, in the simplified 
framework of modelling the electronic repulsion by means of a local nonlinear term. The 
ensuing set of coupled DNSL equations obtained for the vector impurity, is solved numerically 
for arbitrary values of electronic repulsion and magnetic field strength. We focused on the 
cases of a bulk and a surface impurity. In general, we found that it is easier to create 
a bound state in the bulk case than in the surface case. In the latter, the bound state 
diagram showing the number of bound states in terms of nonlinearity and magnetic field 
parameters, show a rather complex structure, that includes a "forbidden region" delimited 
approximately by f/ = B+V and U = B—V. In this regard, the presence of a surface inhibits 
completely the creation of a vector localized mode, in the region f/ ~ 5. In addition, there 
is a parameter region where two bound states exist, although only one of them is (linearly) 
stable. Dynamical results for selftrapping show that, in general, it is easier to selftrap on 
a bulk impurity than on a surface one, although the selftrapping is sharper for the latter. 
This repulsive character of the surface was observed earlier for the one-dimensional nonlinear 



chain 



18|, and seems to be generic to one- dimensional discrete systems. 
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List of Figure Captions 

Figure 1: Left: B — case. Amplitude of localized state for both spins as a function of 
nonlinearity. Right: U — case. Amplitude of localized state for both spins, as a function 
of magnetic field. 

Figure 2: Amplitude of localized mode as a function of nonlinearity and external magnetic 
field, for spin up (left) and spin down (right). Darkest (whitest) shade corresponds to 

Figure 3: (Color online). Spatial mode profiles for U/V — 3 and varying values of B/V. 
Top left: B/V = 1. Top right: B/V = 2. Bottom left: B/V = 4. Bottom right: B/V = 6. 
Solid (dashed) curves denote the spin up (down) mode. 

Figure 4: Left: B — case. Amphtude of localized state for both spins as a function of 
nonlinearity. Right: U = case. Amplitude of localized state for both spins, as a function 
of magnetic field. 

Figure 5: Number of locahzed surface modes, as a function of nonhnearity and magnetic 
field. 

Figure 6: (Color online). Examples of surface locahzed modes. Sohd (dashed) hues refer to 
spin up (down). Top left: B/V = 1, U/V = 4 (stable mode). Top right: B/V = 1, U/V = 4 
(unstable mode). Bottom left: B/V = 3, U/V = 5. Bottom right: B/V = 3, U/V = 1 

Figure 7: (Color online). "Bulk" magnetic impurity: Time- averaged selftrapped fraction 
for both spins remaining on the initial site, as a function of U and B. 

Figure 8: (Color online). Surface magnetic impurity: Time-averaged selftrapped fraction 
for both spins remaining on the initial site, as a function of U and B. 



15 



